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Abstract 

We examine the applicability of the continuum model to describe the surface 
morphology of a hetero-growth system: compositionally-graded, relaxed GeSi 
films on (001) Si substrates. Surface roughness versus lateral dimension was 
analyzed for samples what were grown under different conditions. We find 
that all samples belong to the same growth class, in which the surface rough- 
ness scales linearly with lateral size at small scales and appears to saturate at 
large scales. For length scales ranging from 1 nm to 100 /xm, the scaling be- 
havior can be described by a linear continuum model consisting of a surface 
diffusion term and a Laplacian term. However, in-depth analysis on non- 
universal amplitudes indicates the breaking of up-down symmetry, suggesting 
the presence of non-linear terms in the microscopic model. We argue that 
the leading non-linear term has the form of Ai(V^)^, but its effect on scaling 
exponents will not be evident for length scales less than 1 mm. Therefore, 
the growth dynamics of this system is described by the Kuramoto-Sivashinsky 
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equation, consisting of the two linear terms plus Ai(V/i)^, driven by a Gaus- 
sian noise. We also discuss the negative coefficient in the Laplacian term as 
an instability mechanism responsible for large scale film morphology on the 
final surface. 

PACS numbers: 05.40.+j, 82.20.Mj, 68.55. Jk, 07.79.Lh 
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I. INTRODUCTION 



The dynamics of film growth has proven to be rich and interesting phenomena. In an 
ideal homo-growth, during which the deposition rate is sufficiently low and the temperature 
of the substrate is high enough, the adatoms have enough time to find their optimal positions 
so that most adatoms are registered and the growing front has only small fluctuations around 
the equilibrium shape. The resulting film under this growth condition is smooth. Such 
quality of the film can be maintained indefinitely only if the atoms below the growing front 
are always kept in true equilibrium. In reality, however, growth usually happens in non- 
equilibrium conditions. In fact, the real power of thin film growth is the capability to 
create new materials and to obtain desired physical properties via non-equilibrium growth. 
In practical applications, an often encountered situation is to have the thin film and the 
substrate be different materials, i.e., hetero-growth. In hetero-growth, the growth mode is 
usually not layer-by-layer; p| instead, it depends on equilibrium material properties as well 
as kinetic parameters during growth. In the extreme case, it was demonstrated recently that 
the coherent strain in the film can be utilized to fabricate novel nanostructures. 

During the last decade, much work has been devoted to understanding the non- 
equilibrium film growth. There are several important features observed in the final surface 
morphology. First, surfaces are highly irregular. It is therefore impractical to predict or 
describe such surfaces in microscopic details. A coarse-grained, statistical modeling is more 
appropriate. Secondly, certain large-scale features can survive in the final surface morphol- 
ogy. These features are manifestations of underlying microscopic instabilities. Thirdly, it 
is found, in numerous experiments and computer simulations, that analysis of the surface 
roughness versus the sample's linear dimensions provides a useful classification of growth 
mechanism. Specifically, many surfaces exhibit self-affinity in which a scaling phe- 

nomenon is found 

cr{t) ^ ^{[H{r,t)-{H{r,t))Y) ~ A^L^f [L/m] • (1) 
Here the sample dimension is Lx L, a{t) is the surface roughness, H{r, t) is the height of the 
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surface at position r and time t, and {H{r, t)) is the average height. The length ^(t) denotes 
the characteristic length of the surface. If no other important length scale is present, C,{t) is 
the correlation length built up during the course of film growth and scales as ^(t) ~ (i^t)^ . 
The exponents x and z are useful in characterizing the surface morphology. Different growth 
mechanisms result in different exponents, while the details of growth manifest themselves 
only through non-universal amplitudes A„ and 9. Along with the discovery of the above 
scaling phenomenon, much theoretical work has been devoted to constructing appropriate 
continuum models for describing the film growth. The goal is to reproduce these scaling 
results. 

Experimentally, until now, most work has concentrated on the measurement of the expo- 
nents, which leaves open many important issues. For example, in previous work, no differen- 
tiation in material nature between film and substrate and no examination of the applicability 
of continuum models to surfaces with large features were ever carefully made. Many results 
are thus a priori only applicable to homo-growth systems and to growth without instability. 
Results from model systems |Q show that the quantitative scaling characteristics of surface 
roughness depend on the particular material system; for example, the data for x range from 
0.2 to 1. It is therefore particularly important to examine the applicability of theoretical 
ideas, such as continuum modeling, beyond the homo-growth system. 

In this paper, we examine a hetero-growth system: compositionally-graded, relaxed GeSi 
films on (001) Si substrates. A distinct feature of surfaces in this system is the existence 
of large-scale patterns, known as cross-hatches. These patterns are closely related the 
underlying misfit dislocation network. Much work has been devoted to understanding the 
mechanism responsible for the cross-hatch formation. Here we approach this problem 
differently by examining how surface roughness depends on lateral size, i.e., by a scaling anal- 
yses. In addition, we go beyond the usual approach, that is based solely on the measurement 
of roughness exponents, and perform more comprehensive analyses on non-universal ampli- 
tudes and the up-down symmetry of surfaces. Our results indicate that, up to a length scale 
of 100 /im, a continuum model in which the linear parts are composed of a surface diffusion 
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term and a Laplacian term is appropriate for describing these surfaces. The breaking of 
up-down symmetry shows that non-hnearities must also be present. Detailed analysis sets a 
lower bound of ~ 1 mm for observing the scaling exponents that arise from the lowest-order 
non-linear term. The resulting continuum model that is consistent with the experimental 
data is a two-dimensional Kuramoto-Sivashinsky equation driven by a Gaussian noise. ^ 
We also discuss the instability that might be responsible for the cross-hatch formation in 
the framework of a continuum description. 

This paper is organized as follows. Section II briefly reviews relevant theoretical ideas 
and results. In Section III, we apply scaling analyses to study the surface morphology 
of relaxed GeSi films grown on Si substrates. Both universal scaling exponents and non- 
universal amplitudes are analyzed. The restrictions placed on the proposed continuum model 
by the experimental results and the crossover from this model to other models are discussed 
in Section IV. The Appendixes are devoted to more technical details. In Appendix A, we 
calculate the surface roughness for a sinusoidal surface. In Appendix B, we apply the Dyson- 
Wyld renormalized perturbation theory to analyze one of the possible growth models, that 
can account for our results. 

II. CONTINUUM MODELS AND THEORETICAL RESULTS 

We begin with a brief review of the theoretical situations. There are two useful asymp- 
totic behaviors for f{y) in Eq.(0). When L <^ ^{t), the whole sample evolves "coherently" 
in the sense that a{t) is independent of time; hence f{y) ~ 1 as y ^ 0. It implies that, in 
the limit of large t (L ^ ^(t)), a ~ A^^L^. At the opposite limit, when L ^ ^(t), the local 
surface roughness has not detected the existence of the boundary of the sample, so a(t) does 
not depend on L. One then deduces that f{y) ~ as y ^ oo. That is, in the limit of 
small t , a A„{utf = Att^ , with (3 = xjz. 

In order to explain the above scaling phenomena and calculate the relevant exponents, 
a number of continuum models have been proposed. |TD|-|T^ In these models, H{r,t) is 
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coarse-grained. The interface growth is modeled by noise-driven dynamics 

^ = F[h]+rj{v,t). (2) 

Here h = H — {H), {H) being the average height and 7]{r,t), representing fluctuations in 
the deposition flux, is modeled by a Gaussian white noise with the two-point correlation 
function given by 

(r/(ri,ti)r/(r2, ^2)) = 2r/o5(ri - r2)5(ti - ta). (3) 

The rest of the growth dynamics is lumped into F[h] in which different driving forces for 
adatom movement are represented by different terms. A characteristic of these models is 
that the resulting surfaces are usually driven towards being self-affine, resulting in the above 
scaling behaviors. 

The continuum models can be classified into two classes: conservative and non- 
conservative. In non-conservative dynamics, the flow of atoms onto the surface is as- 
sumed to be normal to the surface. The Kardar-Parisi-Zhang (KPZ) model, in which 
F[h] = u'V'^h + Ai(V/i)^ with z/ > 0, represents the lowest-order realization of such dy- 
namics. Here Ai(V/i)^ accounts the fact that the growth is normal to the interface, and the 



desorption, accounted by uV'^h, is assumed to be important. |12| , [T3| The exponents of this 
model for dimensionality (d) of 1 + 1 are known exactly: x = 1/2 and z = 3/2, while only 
numerical results are available for d = 2 + 1. [|l^] 

In conservative dynamics, one only includes the flow of atoms parallel to the surface. 
Therefore, F[h] must take the form 

F[h] = -V- J, 

where J is the surface adatom current and the total volume / / dxdyh{x, y) is conserved. It 
has been argued that conservative dynamics is the main scenario that occurs in the molecular 
beam epitaxy (MBE) growth. In particular, surface diffusion, rather than desorption, is the 
dominating factor. [|10| -[T2[| The surface diffusion is assumed to be driven by the energetics 
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of adatoms on surfaces so that J obeys the Pick's law J = — aV/i, where /i is the chemical 
potential on surfaces. A couple mechanisms contributing to /i have been proposed. First, 
because adatoms have less (more) bonding opportunities when they reside on "mountains" 
("valleys"), /i is proportional to the curvature. Secondly, adatoms on slopes either have 
higher kinetic energy or their chemical bonds are more stretched. One thus expects that 



they have higher /i. This effect may be accounted by the term (V/i)^. ||TT1 Combining 
these two mechanisms, Lai and Das Sarma (LS) investigated the model in which F[h] = 
— i^V^/i + A2V^(V/i)^. [0 Here —KV^h is the lowest-order contribution from the curvature. 
They found the exponents to be x = 2/3 and 2; = 10/3 for = 2 + 1. 

The existence of a Laplacian term is an important question, for in order to see the 
exponents predicted by the LS model, there must be no Laplacian term. A non-vanishing 
Laplacian term, no matter how small, will make A2 irrelevant, resulting in different scaling 



exponents from the LS results. [|T6| There are several mechanisms that can generate Laplacian 
terms. The typical way to generate the Laplacian term is via the desorption, which is small 
but is not identically zero. [|T| If surfaces are liquid-like, a Laplacian term can arise from 
non-vanishing surface tension. More recently. Villain |T^ argued that the combination of 
surface diffusion on terraces and step-edge (Schwoebel) barriers results in a Laplacian term. 
The coefficient u of the Laplacian term, however, is negative when growing on singular (i.e., 
high-symmetry) surfaces. 

Theoretically, a negative z/ leads to unstable surface growth, producing large-scale fea- 
tures on surfaces. This phenomenon was observed in a recent Monte Carlo study. For 



continuum models, when the negative Laplacian term is combined with the surface diffusion 
term -KV^h and the KPZ non-linear term, the resulting dynamics (without noise), known 
as the Kuramoto-Sivashinsky (KS) equation, is particularly interesting. |^ The KS equa- 
tion has been used to describe a wide class of phenomena associated with instabilities such 
as flame propagation and chemical turbulence. It is linearly unstable, but is non-linearly 
chaotic, exhibiting spatio-temporal chaos. Many years ago, Yakhot ||TB[ conjectured that 
the non-linear term in the KS equation (with random initial conditions) will "renormalize" 



the negative v so that the effective large-scale behavior of the KS equation can be described 
by the KPZ equation with a positive, effective v. In one dimension, this is confirmed by a 
recent work of Chow and Hwa. [|19| In two dimensions, however, conflicting results have been 
reported. There appears to be no consensus about the scaling behavior of the KS 

equation at large scales, though they all agree that an effective, positive v must be present 
at large length scales. Experimentally, quantitative evidence of a non-vanishing Laplacian 
term was not established until recently, although an earlier study used a negative v to 
explain the large-scale features in homo-growth of GaAs films. 

A negative v is not the only possible mechanism for generating large features on surfaces, 
i.e. rough surfaces. In fact, it has been recognized that a flat surface under stress may 
become unstable or metastable. Several stress-induced instability mechanisms have been 
proposed. In the continuum elastic theory, ||2^ an elastic energy is added to the chemical 
potential, resulting in a linear term, cA;^/i(k, t), in Fourier space. Here c is a positive constant. 
Tersoff et al. have proposed a long-range attractive interaction between steps on vicinal 
strained layers. The interaction leads to step-bunching instability, and effectively, it also 
introduces the term /c'^/i(k, t) but with a different coefficient c'. It is important to note that 
these theories were derived for strained films without dislocations. The Ic" term is nonlocal 
and not analytic in real space, and thus is absent in the gradient expansion of F [/;,]. 



III. EXPERIMENTS AND RESULTS 

GeSi films studied here were grown either by MBE or by low pressure chemical vapor 
deposition (CVD). [^] In order to minimize threading dislocation density, the Ge concen- 
trations in these films were increased linearly with the average film thickness until desired 
compositions were achieved. The growth temperature for these samples was sufficiently high 
so as to achieve strain relaxation during growth. The average lattice constants of the films 
are the same as those of bulk crystals, i.e., completely relaxed. Detailed characterization of 
these samples is given in Refs. 6 and 26. We concentrate on samples with approximately 
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the same final Ge composition but grown under different conditions. We also examine sam- 
ples grown with different grading rates. Relevant growth information and non-universal 
amplitudes obtained from our analysis for five samples are summarized in Table I. 

The surface roughness was measured using a scanning force microscope (SFM). Large- 
scale (48 fimy SFM images of all five samples are shown in Figs. 1 (a)-(e). All samples 
except sample C display a long-range ordered cross-hatch pattern on the surface. We 



shall see that the ordering of the cross-hatch pattern does not affect the scaling behavior of 
the surface roughness. In Figs. 2(a) and 2(b), we show samples B and C respectively, on a 
magnified scale, (14.5 /xm)^. The line cuts (height changes versus lateral distances) indicated 
on Figs. 2(a) and 2(b) are shown in Figs. 2(c) and 2(d), respectively. To determine the 
surface roughness at length scales varying from 1 nm to 100 fim, images (256 pixel x 256 
pixel) similar to Figs. 1 (a)-(e) were taken at different length scales and on at least two 
random spots for a given sample. When performing scaling analysis, we divided each image 
into smaller images ( 128 x 128, 64 x 64, 2 x 2). The roughness for a given length scale L 
was obtained by first calculating the rms roughness inside each L x L image, then averaging 



over the ensemble of images of the same size. |23] To obtain data for L from 1 nm to 100 /im. 



we combine results from several images of different sizes. 

Fig. 3(a) shows cr versus L in a log-log plot for all five samples. These curves all show 
similar behaviors, i.e., a oc at small length scales over 3 decades in both a and L, and 
a flattens out with zero slopes {x ~ 0) above some sample-dependent crossover length 
Lj, = 0.5 ~ 10/im. At larger length scales, all surfaces are fiat in the sense that the surface 
roughness reaches a saturated value (cxo) that is almost independent of the system size. For 
the five samples we studied, ctq ranges from 37 to 186 A, and L,y from 0.7 to 5.5 jjm (see 
Table I). These differences are more easily seen in Figs. 2(a)-(d), especially the difference 
for (To [Figs. 2(c) and 2(d)]. Note that the y-axis (height change) in Fig. 2(d) is twice 
that in Fig. 2(c). When we re-scale L by L,y and a by ao, all the data collapse onto one 
universal curve, as shown in Fig. 3(b). The successful collapsing demonstrates that, despite 
the differences in synthesis parameters and apparent surface morphology, all samples belong 
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to the same growth universahty class governed by (t/uo = f{L/L^), where the collapsed 
data points in Fig. 3(b) trace out the universal scaling function /. The collapsing of data 
indicates that all samples can be described by the same continuum model, if there exists 
an appropriate one. It also implies that the order of the cross-hatch pattern is not relevant 
since the surface of Sample C does not display the ordered the cross-hatch pattern but can 
still be collapsed with the other data. 

A simple power counting suggests that scaling results shown in Fig. 3 can be accounted 
by the linear model: 

F[h,t] = uV^h-KV^h. (4) 
Since the —KV^h term represents surface diffusion []T2|,pP[ which tends to stabilize growth, 



K is taken to be positive. In this model {d = 2 + 1), when —K'V'^h dominates, x = 1; while 
if u'V'^h dominates, x = 0. |jl| The crossover length L^, between the two scaling regimes is 
setby27r(fs:/|z/|)i/2_ ^he good agreement in scaling exponents between experimental results 
and the model shows that surface diffusion dominates at small length scales. 

In a more detailed analysis, we introduce an infrared cutoff ko{= 27[/L) and an ultraviolet 
cutoff A (= 271 /a, where a is the lattice constant) so that the roughness can be calculated 
via the integral 

^(L) = £ / (Mk, ^) M-k, - uj)). (5) 



When h dominates, we find o-j, ~ ^ rjo / {2tt\i'\) y\n{L/a). In Fig. 4, we show our fit for 
sample B. The lattice constant read off from this fitting is about 6A, which is comparable to 
the lattice constant of Geo-sSio.?. Similar fittings have also been done for the other samples, 
and they are all consistent with this form. In this regime, the non-universal amplitude A^^^ is 



r7o/(27r|z/|). The experimentally observed values for A^^'^ ranges from 12 to 65 A(see Table 



I). The saturation roughness ctq is related to A'^^^ via the relation: ctq = A^^\/\n{Lo/a), where 



Lq is the scan size in the Laplacian term dominant region. Since for all samplesYln(Lo/a) is 
about 3, ctq and A^^ are of the same order. Using L^, determined from the crossover of the 
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scaling exponents, and A^J:'\ determined from a logarithmic fit of the saturation region, to 



calculate A^Jf^ Jlii{L^/a), the results agree very well with experimentally measured cxo for all 



five samples (see Table I). The reason A^^'' ^\n{Lu/a) is consistently smaller than o"o comes 
from the fact that L^, < Lq. At small length scales (L < L^y), if —K'W^h dominates, we get 



<^K ~ y'r]o/{^Q'^^K)L. Here A^^^ is ^Jrjo/{lQTi^K), which is the slope da/dL for L < L„ . 
Experimentally, we can independently determine 3 parameters: (Jq, A^/"^ and Lj,. For our 
data to be consistent with Eq.(|), these parameters must satisfy ao/A^J^^ ^ A.75Ly. In Fig. 
5, we plot ctq/A^^^ versus L^. We see that the linearity between (Jq/A'^J^^ and is quite 
good, though the slope is 1 and not 4.75. This shows that as far as scaling is concerned, 
the above linear model describes the growth of this real, complex hetero-growth system 
reasonably well. 

The discrepancy between the measured slope for ao/A^J^'' versus L„ and the predicted 
value from scaling suggests that non-linear terms do not identically vanish. Therefore, as a 
test of linearity, we also examine the value of = \ {h^)\^^^. In all samples, we found nonzero 
6s (see Table I). There is no apparent enhancement of 6 from increasing the grading rate 
(samples B and C). In Fig. 6, we show the plot of 6 versus L for sample B. Even though the 
data are more noisy, it shows that 6 has a similar dependence on L as that of the roughness. 
This indicates that these surfaces are not multi-affine. The non-vanishing of 6 implies 
that the up-down symmetry is broken and thus a non-linear term must be present. The 
lowest-order non-linear term appears to have no effect on the scaling behavior of a below 
100 /xm (including the linear relation between ao/A^J^^ and L,y), but they are important in 
determining factors besides scalings such as the slope between ao/A^J^^ and L^. 

We can further estimate the magnitudes of u, K, and rjQ. Assuming that surface diffusion 
is induced by variations of the chemical potential, K is given by Dg'-fQ'^n/kBT, p9[ where Dg 
is the surface diffusion constant, n is the areal density of adatoms, 7 is the surface tension 
and Q is atomic volume. At relevant growth temperatures, Dg is of order 10~^cm^/sec, 
3T| while 7 ~ 10^ erg/cm^ . The values of the remaining parameters are standard. ||32 
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Altogether, we find that K ~ 10~^'^cm^sec~^ , which implies that \v\ is of order 10~^/im^sec~^ 
via 27r(_ft'/|z/|)^/^ ~ Lj, ~ lyum. From the relation A^^^ ~ ^J rj^ / {2tx\i'\) and the experimental 
values of A^^'^ and |z/|, we find r^o ~ 10^^^ cm^/ sec. 

Another important feature in Fig. 3 is that samples with faster grading rate (80%Ge//im), 
such as samples B and C, are rougher (with larger a and A^J^'^) for L < L^. In addition, a 
for samples A, D, and E, all of which were grown at 10%Ge//im grading rate, are identical 
at small length scales. For a relaxed film, strain fields are not uniform in the films, but 
concentrate near dislocations. |^ The larger grading rate means that the growth surface is 
closer to the dislocations and therefore the surface strain fields are larger. Hence, the larger 
A^^^s observed in the 80%Ge//im grading rate samples suggest that strain enhances surface 
roughness. From a different point of view, since the primary effect of the —KV'^h term is 
to decrease roughness, our results indicate that strain fields suppresses adatom diffusion on 
the surface. 

IV. DISCUSSION AND CONCLUSIONS 

We now examine our results more closely. First, the relation A'^^^A^J^^ ^ has a natural 
geometrical meaning. It reflects the wavy nature of the surface morphology in our samples: 
if one treats A^^^ approximately as the amplitude of the wave, A^J^'^ as the slope from the 
valley to the peak, and Lj, as the width from peak to peak, the relation A^^yA^'^ ^ 
follows from the definition of slope. In Appendix A, we calculate the surface roughness of 
a continuum sinusoidal surface. The surface roughness of this model surface has a similar 
saturation for the L ^ L^, but in the opposite limit L <^ L^, the roughness does not have 
exact linear scaling relation with L. One should be further cautious that the real surface is 
not exactly a sinusoidal wave. In Figs. 2(c) and 2(d), we show the surface roughness line 
cuts for samples B and C, respectively. The local slopes vary from < 10 A/jjm to > 800 
A/yum for sample B and > 1000 A//im for sample C, indicating that these surfaces contain 
more than one wavelength. 
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Secondly, the experimentally observed v (= ve) is an effective one and must be positive. 
If V is negative, /i(k, t) grows exponentially. The roughness would not saturate at large 
times. In this case, one finds the roughness o"(t, L) by expressing the height in terms of r/ 



/i(k,t)= rrft'exp(|z/|A;¥)r7(k,t') 



and calculates the correlation function (/i(k, t)/i(k', t)). The roughness is then determined 
by the integral / rf^k / d2k'(/i(k, t)/i(k', t)). We obtain 



27r|z/| 



2V« exp(2|z/|Pt) , L 
; dk — In — 

2tt/l k a 

EifSlz/k^lVEi (8W\n''l:]-2\n- 



) V L'^ J a 



(6) 



where Ei is the exponential integral function and we have assumed that initially the surface 
is fiat, i.e., cr'^{L, 0) = 0. For our samples, the growth time is about 10'' sec and A ~ lA ^ ^ 
10^(27r/L), so the order of 2\ue\AH is 10^°, which is very large. As a result, the first term 
dominates so that a{L,t) is approximately independent of L. Therefore, if ue is negative, 
its value, when combined with Eq.(^ and the observed ctq, yields a value of r/o at the order 
of 10^^°^cm^sec~^, which is unreasonably small. It therefore implies that ue can not be 
negative. 

What is the bare growth equation that can generate the observed cross-hatch pattern 
and yet result in a positive ue at large length scales? At first, it seems that the term 
ck^hCkjt), induced by the stress, must be present. However, as mentioned earlier, this term 
was derived only for strained films without dislocations, while our samples are completely 
relaxed. Therefore, we do not include it and consider the case when the bare u (= uq) 



that appears in the growth equation is negative. ||12[ This coefficient z/q enters the roughness 
cr{L, t) at very early times. When t is very small, because the height h of the surface is still 
small, the non-linear terms, which are higher order in h, can be neglected. Therefore, F[h] 
in the growth equation can be simply approximated by a negative Laplacian term z/qV^/i. 
Eq.(B|), with u replaced by z/q, then describes the roughness only when t is very small. The 
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value of z/q may be obtained, for example, by measuring the first and second time derivatives 
of surface roughness via the relation 

= |z.o|[A^ + (27r/L)2] ^a\L 

where we have expressed r^o in terms of the initial increase rate of a'^{L,t). After a short 
period of initial growth, the height h has grown so large that one cannot ignore the effects 
of non-linearity any more. According to the standard picture of pattern formation, |^ non- 



(7) 

t=0 



linear terms then saturate the initial unstable growth, resulting in the final morphology. If 



z/q is negative due to Schwoebel barriers, [|12[ we expect the growth on vicinal surfaces will be 
much smoother because the bare u is positive even at early times. This is indeed observed 
experimentally. 

There are a couple of important issues that need to be addressed: (1) What are the 
large-scale |Q scaling behaviors of the final morphology? and (2) what are the non-linear 



terms that enter the growth equation? From our results, it is clear that up to lOO/im, the 
scaling behaviors at large length scales are captured by the growth equation composed of a 
Laplacian term with an effective, positive v. This fact also provides us some insight about 
the leading order non-linear term. In Appendix B, we examine in detail the combined 
growth equation, F[h, t] = uoV'^h — KV^h+ A2V^(V/i)^, within the framework of Dyson- 
Wyld renormalized perturbation theory. |j20| We show that uq does not get any correction 
that comes from A2V^(V/i)^ in the large length scale limit. This results reflects what we 
mentioned earlier: A2 is irrelevant in the presence of a Laplacian term. Thus, we are left 
with Ai(V/i)^ as the only possible leading non-linear term. The resulting growth equation 
is the KS equation, driven by a Gaussian noise. As discussed in Section II, the situation 



of theoretical work on the large length scale scaling behaviors of the KS equation P, p!8|JT9 
is not clear now. In fact, the reported results for the two-dimensional case are conflicting. 
Nevertheless, they all agree that an effective and positive u must present at large length 
scales. This appears to be precisely what we observe. Although we are not able to resolve 
the theoretical conflict using these experimental results, detailed analyses of our results 
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set a minimum size of the sample for resolving this conflict. Because we did not observe 
the scaling exponents predicted by the KPZ equation, if Ai(V/i)^ survives in the effective 
equation, its effect must be small for length scale below 100/im. The lower end of this non- 
linear term dominant regime must match with the upper end of the regime dominated by the 
Laplacian term. Since the characteristic timescale in the non- linear regime is z/L^, matching 
the timescales in two regimes defines a crossover length Lc 

/^N 1/(2-^) 

, (8) 

where z and v are the corresponding exponent and non-universal amplitude in the KPZ 
regime. Another relation can be obtained by matching the roughness. We find that 



^ AHln(Wa) 



where A^^ and At (in the following equation) represent the corresponding non-universal 
amplitudes in the KPZ regime. Since At = A„{i))^, we obtain 



Lc> 



Ai''hn{L,/a)/AtY^\ (10) 



Because we do not observe any crossover below lOO/xm, has to be larger than lOO/im. In 
addition, by setting Lc on the right hand side of Eq.([lO|), we can get an estimate of the lower 
bound for Lc- Using the following figures: z/x ~ 4, A^*^^ ln(Lo/a) ~ lO^A for Lq ^ lOO/im, 
At ~ lA/(sec)^ and u ~ 10~^/im^sec~^, we find that the lower bound is about 1mm, which 
is beyond our measurement range. 

In conclusion, we have observed a universal scaling behavior for the surface morphology 
of compositionally-graded, relaxed GeSi/Si(001) films. Quantitative analyses on scaling 
exponents and non-universal amplitudes show that the scaling behaviors for samples grown 
under different conditions all belong to the same universality class, which can be described 
by the linear model F[h,t\ = u'V'^h — K'V'^h for Inm ^ L ^ 100/i m. In combination 
with further theoretical analyses, it is argued that the underlying growth model is the KS 
equation driven by a Gaussian noise. Our results indicate that as far as the roughening 
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exponents are concerned, up to 100 //m, the effective tfieory of the KS equation is the hnear 
model, pushing the length scale for observing the KPZ scaling exponents to be above 1 mm. 
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APPENDIX A: 



In this appendix, we demonstrate some features of the surface roughness for a continuum 
wavy surface. Specifically, we shall consider a sinusoidal surface y — yosin(27rx/Lo), but 
many features do not depend on the particular form we choose. Since one only measures 
discrete points in experiments (the number of pixels determining the distance between adja- 
cent points), we assume that Lq — NqE, where A^o is a positive integer and e is the distance 
between adjacent points. The experimental data then consist 2A"-|-1 equidistant points (a;„, 
hn ) centered at {xo, ho), where we have defined hi = yosm{27:Xi/ Lq) and Xn — xo = ns with 
n ranging from —N to N. The surface roughness can be found by evaluating the following 
sums 



{h) 



yo 



2N + 



2N + 



n=N 


o 




n=-N 




n=N 




E 


sin 


n=-N 





27r(n£ + Xo) 
Nos ' 

27: {ne + Xq) 



(Al) 
(A2) 



For large Nq, these sums become integrals. Evaluating these integrals, one obtains a surface 
roughness a{N, Lq, xq), which depends on both xo and Lq. If the modulation on the surface 
contains only one wavelength, the global surface roughness is simply an average of xq over 
the period Lq. We find that 



yo 



L 



f Lc 



1 — cos ■ 



AttL 



nl/2 



Lq J 



(A3) 



2L + e V2L + £27ry 

where L {=Ne) is the size of the sample. In this case, a approaches yo as L approaches oo. 
For L < Lq, a does not scale; instead, it oscillates with L. 



APPENDIX B: 

In this appendix, we shall investigate the large scale behavior of the model 
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8h 

— = uoVh - KV^h + A2V' {Vhf + r/, 



(Bl) 



where is negative. Specifically, we shall show that there is no correction to z/q at large 
length scales. 

Two important functions we shall work with are the response function G{k.,uj) and the 
two-point correlation function U(/c, a;). They are defined by 



r^d. ._/Sh{'k,uj 



\J{k,uj) = {h(k,uj)h{-k,-uj)) . 



(B2) 
(B3) 



The generic forms of the these function are 



G{k, uj) = — — ■ — — , 

\J{k,ij) = \G{k,u)m27]o + <l>{k,u)] 



(B4) 
(B5) 



where 7^ = —v^k"^ — Kk^., J]{k,uj) and ^{k,uj) are the self-energies. We shall be interested 
in scaling solutions in which G{k,uj) = g{uj / uk^) /vk^ and \]{k,uj) = u{ijj/uk^)/k'^. These 
scaling solutions are asymptotically correct only for A; ^ and u; ^ 0. Therefore, 7^ has 
to be subdominant to vk^ for A; — 0; hence z <2. The exponents z and A are related to x 
via the relation: 2% = A — d — 2;. We shall show that z = 2 is not acceptable. 

Let us first consider the one-loop diagrams in the renormalized Wyld-Dyson perturbation 
theory. |20| The self-energy is 



q Jn 



A;2(k-q)-q g^k ■ (k - q) G(g, l])U(|k - q|, - fi). (B6) 



By appropriate changes of variables: Q = q/fc, t = Vt/vq'^ , and s = u/uk^, we can rewrite 
Eq. (|B6|) in the following forms 



S(2) {k, u) = k''+^'^ a(') {s, A/k, mo/k) . 



(B7) 



Here a^"^^ is a function of s, A/k, and mo/k. A and mo are the ultra-violet (UV) and 
the infrared (IR) cutoffs in the g-integrals. These forms imply that if the g-integrals are 



divergent, the leading terms in Ti^'^\k,uj) must take the form ^ ^A*^ or k'''^^ 



m, 



; 



depending on whether the divergences are UV or IR. Here both S and 6' must be positive. 
If the divergences are UV type, we set all internal momentum, such as q and |k — q|, to 
A, and all internal frequencies to A^. Only external momentum are left intact. Since in 
S*^^^(A;, uj) the first vertex carries one k'^ and the final vertex carries one external momenta k, 
the lowest term would be 0(A;^k). However, because S(fc, cj) depends only on k, this term 
must vanish. Hence, the leading term is k^A^~^'^~^. The sub- leading terms are k^A'^~^'^~^, 
^8^d-A ^ g^^^ Qj^^ with the power of A decreasing until the power of A becomes negative. 
For the above to be correct, one requires A < d + 4, and the difference between A and 
d + 4 decides the number of sub-leading terms. It is clear that the UV divergences do not 
contribute any correction to Uq. 

On the other hand, if the divergences are IR type, one sets g or |k — q|, but not both, 
to uiq. Simple power counting leads to the conclusion that the leading term of T.^'^\k,uj) 
is k^~^ {moY~^^~'^'^^ . The sub-leading terms are terms with less power of k, so they could 
correct uq. The correction must take the form k'^(moY~^^~^. Obviously, it implies that 
A > d + 6, which results in x > (6 — 2;)/2 > 2. Since in the physical regime, x ^ 1) 
this is also ruled out. Therefore, the IR divergences do not contribute z/q in the physical 
regime(x < 1) either. 

Finally, if the integral cr'-^^ is convergent, it depends only on s in the limits A ^ oo and 
mo 0. In order that Uq gets a correction, we require T.^'^\k, cu) ~ k'^. Hence d + 8 — A = 2, 
i.e., A = d + 6. However, this values falls into the regime where the integral a^'^^ is not 
convergent, but IR divergent! 

The above analysis can be easily generalized to higher order terms. We find that the 
dimension of 2nth order terms in J2 is n{d + 8 — A — z) + z. These terms can not be both 
convergent and at the same time contribute 0{k'^) because it would imply A = d + 8 — 
2/n + {1/n — l)z, which is greater than d + Q and thus falls into the IR regime. Thus these 
terms must be divergent. We find that the leading contribution of the UV divergences to 
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E are of the form OnA"^'^"'"^ ^, where a„ is the contribution of 2nth order terms; 

n=l 

hence the UV divergences do not correct uq at all. Similarly, for IR divergences, possible 
corrections to coming from the 2nth order terms must be of the form k'^rn^^'''^^ ^ ^^^^ ^. 
It implies that A > d + 8 + (1/n - l)z - 2/n, and hence x > 4 + (l/2n -l)z -l/n> 2, 
which is not in the physical regime, so the IR divergences do not contribute vq at all in the 
physical regime(x < 1). We thus conclude that to all orders in the Wyld-Dyson renormalized 
perturbation expansion, there is no correction to vq. 
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TABLES 

TABLE I. Detailed growth parameters {Ri and R2 are the growth and grading rates) and 
non-universal amplitudes obtained from our analysis of samples A - E : ctq and 60 are saturated 
values of surface roughness and , respectively. 



sample 


A 


B 


C 


D 


E 


Ge(%) 


30 


30 


30 


40 


30 


growth method 


MBE 


MBE 


MBE 


CVD 


CVD 


T(° C) 


900 


900 


900 


850 


650 


P(mT) 


- 


- 


- 


2 


50 


i?i(A/scc) 


3 


3 


3 


7 


12 


R2{Ge%/nm) 


10 


80 


80 


10 


10 




1.5 


0.67 


1.05 


2.0 


5.5 


^0 (A) 


37 ±4 


65 ±4 


186 ± 15 


61 ±7 


124 ±5 


^o(A) 


13 ±4 


19 ±5 


91 ±38 


29± 11 


51 ±6 


A^^^ (A) 


12 


22 


65 


20 


38 


Ai""^ (10-4) 


26 


87 


169 


29 


28 


At\/hi(L,./a) (A) 


31 


58 


178 


57 


115 
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FIGURE CAPTIONS 



FIG. 1. SFM images of samples listed in Table I: A(a), B(b), C(c), D(c), and E(d), 
respectively. The image sizes are (48/im) x (48/im). 

FIG. 2. (14.5 //m) X ( 14.5 //m) SFM images of samples B(a) and C(b), respectively. 
The surface roughness line cuts (height change versus lateral distance) indicated by the 
white lines in (a) and (b) are shown in (c) and (d), respectively. 

FIG. 3. (a) Surface roughness (u) versus L for all 5 samples in Table I. The calculation 
of cr(L) is given in the text. The line represents cr cx L^. (b) cr/uo vs. L/ for all samples. 

FIG. 4. A logarithmic fit to the data of sample B. The lattice constant obtained from 
this fit is about 6A, and A^^^ = ^rjo/ {2t:\p\) is 22 A. 

FIG. 5. Ratio of non- universal amphtudes gq/A^^^ vs. the crossover length L^,. The 
line is a guide to the eye with slope = 1. 

FIG. 6. The third moment 9 = versus L for sample B. 
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